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We derive nonlinear o models (chiral Lagrangians) over symmetric spaces U(n), U(2n)/Sp(2n), 
and U(2n)/0(2n) from U(A r ), O(N), and Sp(2iV) lattice gauge theories coupled to n flavors of 
staggered fermions, in the large-iV and g 2 N limit. To this end, we employ Zirnbauer's color-flavor 
transformation. We prove the spatial homogeneity of the vacuum configurations of mesons by 
explicitly solving the large-iV saddle point equations, and thus establish the above patterns of 
spontaneous chiral symmetry breaking without any assumptions. 
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Introduction. The concept of spontaneous symmetry 
breaking introduces an essential parallelism between the 
spectral statistics of disordered condensate Q and the 
low-energy dynamics of strongly coupled gauge field the- 
ories [||, involving a juxtaposition of the lightest and 
weakly coupled Goldstone particles in these theories: dif- 
fusion modes associated with spontaneous breakdown of 
the symmetry between advanced and retarded Green's 
functions, and pions associated with spontaneous break- 
down of the chiral symmetry. However there is a con- 
siderable difference at the practical level. An ensemble 
averaging can be analytically performed for disordered 
Hamiltonians, so that the Goldstone manifold and the 
nonlinear a model over it describing the low-frequency 
regime of the theory manifest themselves without ambi- 
guity. On the other hand, gauge theories with propagat- 
ing gluons do not admit such an analytic computation. 
Accordingly, in order to identify a low-energy effective 
theory one needs to appeal to nonperturbative theorems, 
namely Vafa-Witten's theorem M that forbids sponta- 
neous breaking of the vector symmetry in a vectorial the- 
ory, and 't Hooft's anomaly matching B between funda- 
mental and effective theories. Relying on these theorems, 
non-standard patterns of chiral symmetry breaking have 
been studied || in the context of technicolor models. 

The symmetry-based approach puts a strong restric- 
tion on operators that are allowed in the effective chiral 
Lagrangian, but the setback is that their coupling con- 
stants are left as free parameters that can by no means be 
related to those of the fundamental theory. To overcome 
this difficulty, attempts |(| were made to derive an effec- 
tive theory directly from a microscopic model under an 
extreme condition, namely a U(N) lattice gauge theory 
in the strong coupling limit, combined with the large- N 



limit. This serves as a basis on which a strong coupling 
expansion is performed and the plaquette action can be 
taken into account perturbatively in (g 2 A^) _1 

Recently Altland and Simons || have renewed interest 
in this subject. Instead of computing a one-link integral 
via the Brezin-Gross formula (^] valid only in the large- 
N limit, they have adopted an alternative method, Zirn- 
bauer's color-flavor transformation [ ^0|JTl| ], which has a 
clear advantage of being valid at finite N. By converting 
the integrations over strongly coupled gluonic variables 
into supposedly weakly coupled mesonic variables, they 
have provided a transparent derivation of the low-energy 
effective action corresponding to fermions in the complex 
representation. However, their and all the previous works 
have relied on a crucial and far-from-evident assumption 
that the mesonic fields in the vacuum are homogeneous. 

The purpose of this paper is twofold: by applying 
the color-flavor transformation to strongly coupled lat- 
tice gauge theories with three classical gauge groups and 
with staggered fermions, we shall i) demonstrate the spa- 
tial homogeneity of the mesonic fields in the vacuum, 
thereby establishing the chiral symmetry breaking rigor- 
ously, ii) extend the method of Altland and Simons valid 
for quarks in the complex representation of the gauge 
group, to other two classes of representations, so as to 
derive non-standard type of chiral Lagrangians Jl^| . 

Lattice Gauge Theory with Staggered Fermions. We 
consider a lattice gauge theory coupled to staggered 
fermions on an even(= d) dimensional lattice Z d 3 x, 
which is split into odd and even sublattices labeled by o 
and x , respectively. The set of 2d odd (even) sites adja- 
cent to an even (odd) site x will be denoted as o x (x x ). 
We write the partition function with an emphasis on the 
bipartite nature (the lattice spacing is set to unity): 



Z(g; m) = J dU xy J dip x dip x exp 



S[U] + £ 



fix 



(i) 



i 



Here the link variable Ul* takes its value in the 
gauge group G = U(N),0(N), or Sp(2A). The 
Haar measure of G is denoted by dU, and S[U] — 
9~ 2 T,{ xy zw) tTU *v U zy U ™ U L + c - c - is the _ plaquette 
gauge action. The site variables ip^ a and 4> x a are m ~ 
dependent Grassmannian numbers. The color indices 
... run from 1 to N for G = U(N), 0(N) or to 2N 
for G = Sp(2JV), and the flavor indices a,b, . . . run from 
1 to n. The phase is defined by r\ xx ±f L = (— l) 2 ^ \ 
and the quark mass matrix is by m = diag(mi, . . . , m n ). 
As the fundamental representation of 0(N) (Sp(2A)) 



is real (pseudoreal), i.e. its symmetric (antisymmetric) 
product contains an invariant, the chiral symmetry group 
is extended to U(2n) for these cases §. This can be made 
apparent by introducing 2n-flavored site variables 



i.A 



ty? ft*] or £^ a ] (y G o), (2) 
[$r a V£ a ] or (sex), 



(A, B = l,...,2n,E = «j 2 <g> i w ) for G = 0(N) or 
Sp(2A), resp., and writing the fermionic action as 



E ^^O^-iE^^^ + sE^'^^^ (O xy eO(N)), (3a) 

E iE v """ A " ; '" ,,il "« ; " • ^E , "'" A '"'" Uiv "' ; " (^esp^). (3b) 
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Extended mass matrices m in the above are defined by 
m = («<T2 or (7!) ® diag(mi, . . . , m n ) for G = O(N) or 
Sp(2A), respectively. The protected vector subgroup of 
this extended chiral group U(2n) is read off from Eq. (^) 
to be Sp(2n) for G = 0(A), or 0(2n) for G = Sp(2JV). 

Color-Flavor Transformation. Henceforth we concen- 
trate on the strong coupling limit g 2 N — > oo. The theory 
is thus referred to as a non-Abelian random flux model. 



In this case, r\ xy can be absorbed into the redefinition 
of the link variables. We can reexpress the partition 
function of this model in terms of an integration over 
flavor-singlet link variables into the one over color-singlet 
link variables, by Zirnbauer's color-flavor transformation. 
That for the U(A) link variable coupled to fermions reads 
up to an irrelevant numerical factor, 



dU exp 

U(AT) 



dZ ab dZ ab 



rii<a,&<i 

det Ar + 2n (l + ZtZ) 



exp 



U\' a Z ab i/;i' b + ^j°Zt°fyo 6 ) . (4s 



The integration on the RHS is over complex n X n matrices, Af = GL(n, C). Although the case with the O(JV) link 
variable coupled to fermions was not explicitly provided (see Ref. Jllj for its bosonic counterpart), one can show that 



/ dO exp(¥' A O i:i ^ j ' A ) = f 

Jo(N) JZ 



III 



<o<b<2n 



dZ ab dZ ab 



z= _ z t det^+^-^l + ZtZ) 



exp (&< A Z AB ¥< B - W> A Z^ AB W> B ) 



(4b) 



The integration on the RHS is over the set of complex 2n x 2n antisymmetric matrices, which is diffeomorphic to 
J\f = GL(2n, C)/Sp(2n, C). The case with the Sp(2iV) link variable coupled to fermions reads 10 



dS exp Vs'mI' " M 

Sp(2JV) 



ni< 



z=z t det N + 2n + l (I + ZiZ) y \2 2 1 



The integration on the RHS is over the set of complex 2n x 2n symmetric matrices, which is diffeomorphic to 
Af = GL(2n, C)/0(2n, C). Applying these transformations to each link and then integrating over fermions, we obtain 



z(m) = i n 



dZ X y dZ X y 



N iex, yGo det 



(1 + Z xy Z x y) x 



n det^ E^- 2m n detN e z i * 



(5) 



Here e — 0, —1, 1 for G = XJ(N), O(N), Sp(2A), respectively, and we have redefined N/2 — > N and m — > 2m for the 
0(A) case. Pfaffians from fermionic integrations Pf(A) = J d^ exp(^E'A4') arising in the 0(A) and Sp(2A) cases 
are replaced by det X / 2 {X). In spite of residing on a link, Z xy variables resemble the mass matrix that couples to a 
fermions bilinear on an identical site, thereby indicating the spontaneous chiral symmetry breaking. 

Homogeneity of the Saddle Point. So far the manipulation applied to our non-Abelian random flux models was 
exact. We finally make an approximation, by taking the N — > oo limit. The saddle point equations read 
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Z *y = E Z l'v ~ 2m > ^Ur 1 = E Z *V ~ 2m ( 6 ) 

a;' g X y y' £oj; 

(x G x, y G o x ) for G = U(iV). Those for G = 0(JV) or 
Sp(2A r ) are identical as the above, with Z xy constrained 
to be antisymmetric or symmetric. 

Now we shall demonstrate that the unique solution to 
Eq.(Q) at m — is given by 

Z xy = {2d-l)- 1/2 U, U G U(n) for y x G x, v y G o x , (7) 

i.e. the vacuum configuration of the mesonic variable is 
homogeneous. In other two cases, we simply need to 
restrict U to be an antisymmetric (G = O(iV)) or a sym- 
metric (G = Sp(2iV)) unitary matrix, which in turn can 
be parametrized by another unitary matrix V G U(2n) 
by the canonical projection U — VJV T (J = ia 2 £g> 1„) 
or U = VV T , respectively. 

The set of saddle point equations involving an even 
site x read, after using abbreviations Z p = Z XVv for 
2/i, •• ■ ,U2d G o x , 

2d 

(K)~ 1= E ^ P=l, ...,2d(>4). (8) 

9=1 («#p) 

Here Z 1} . . . , Z 2d G GL(n, C). Eq.(|) clearly leads to 

(Zl)- 1 + Z p = (Zl)- 1 + Z q , Vp,Vg = l,...,2d. (9) 

We make singular value decomposition of the matrices 

Z P = U P S P V P , S p = diag(s ( f ) ,...,s^), (10) 

where U P ,V P £ U(n) and , . . . , si^ > 0, and define 
ci p ' =1/ si P ' + ■ Without loss of generality we choose 
(U p , V p ) so that <t[ p) >■■■> cri p) . From Eq.(|), we have 

u p (Sp 1 + s p )Vp = u^s- 1 + S q )V q 
=> up^s; 1 + Sp)VpV} = v q v p \s p l + S p )U}U q 

= S- 1 + S q (11) 
=> UlU p (S- x + S p fulU q = (S- 1 + S q f. (12) 

From the uniqueness of the eigenvalues, we find that the 
diagonal matrix Sp -1 -I- S p is independent of p, and thus 
is simply denoted as Sp -1 + S p = diag (a\, . . . ,a n ) . Then 
Eq.(|l2"|) in a componentwise notation reads 

(UlUp) ab (a 2 b - o*) = 0, 0,6=1,...,!*. (13) 

As Cb + a a is nonzero, we find 

(Up p ) ab (a h -a a ) = 0, o,6 = l, ...,n (14) 
O U^UpiS- 1 + Sp) = (Sp 1 + Sp)UlU p . (15) 

Using Eq.(la) in Eq.(Q), we obtain 



Upv = V o V p = W ip- ( 16 ) 
Note that W^ p = W pq . In terms of W qp , Eq.(||) reads 

2d 

S p 1 = E W 1 S i W «pi P=l,...,2d. (17) 

9=1 (9#P) 

Let us now suppose <7i = • • • = oy > oy+i > • • • > a n 
and consider possible two cases separately. 
In the case r = n, we have 

Sp 1 + S p = al n , p=\,...,2d. (18) 

We rewrite the saddle point equations (0) as 

2d 

WjpS r W rp = S; 1 - W^S^p, (19) 

r— 1 (r^p.q) 

where p, g = 1, . . . , 2d and p ^ q. Since W? p S r W rp are 
positive-definite, the LHS of Eq.(|l9|) is positive-definite. 
Therefore S*" 1 and Wl p S q W qp are ordered hermitian ma- 
trices []l3| so that the largest and smallest diagonal ele- 
ments (eigenvalues) of Sp 1 are larger than the largest and 
smallest diagonal elements of S q , respectively. Suppose 
the largest diagonal element of Sp" 1 is s. From Eq. (fl8|) , 
the diagonal elements of S p and S q are s or 1/s. The 
largest diagonal element of is 1/s and s > 1, since it 
must be smaller than s. Therefore all the diagonal ele- 
ments of S q are 1/s. Then the smallest diagonal element 
of S^ 1 has to be s, because it must be larger than the 
smallest diagonal element of S q . Consequently all the 
diagonal elements of Sp 1 are s and we conclude 

Z p = UpSpVp = UpVpSp = U q V q S q = U q S q V q = Z q . (20) 

In the case r < n, we obtain from Eqs.Q and ©, 

cw^)- = o far { « z i'+ i;. r :. , fc „r iV:::" (21) 

for an arbitrary pair of indices p and It means that 
the saddle point equation (0) is decomposed into those 
with smaller ranks r and n — r. Thus we can inductively 
show Z\ = ■ ■ ■ = Z 2 d, using the argument given in the 
previous case for each irreducible component. 

Substituting Z\ = ■ ■ ■ = Z 2 d into Eq.(^), we obtain 

Z p = (2d- iy 1/2 U, UeV(n), y P =l,...,2d. (22) 

As this proof applies for each of even sites, or odd sites 
with Z xy substituted by Z\ y , we finally establish Eq.(^|). 

Chiral Lagrangian. Next we take into account a small 
deviation from the chiral limit, as well as fluctuations of 
the mesonic field around the vacuum configuration that 
has been proven in the last section to be homogeneous. 

The terms to be collected are of lowest nontrivial orders 
in masses and momenta, which are 0(m 1 ) and 0(d 2 ), 



3 



respectively. In the case of G — U(iV) B, we split the 
mesonic variable into massive and Goldstone modes as 



Z xy = (2d-l) 



-l/2p IT 



(23a) 



where P^j, € GL(n, C)/U(n) and C/^j, € U(n). In the 
case of G — 0(A) or Sp(2A), the mesonic variable is 
constrained to be antisymmetric or symmetric, respec- 
tively, so we must employ parameterizations 

Z xy = (2d - L) ; 2 l>r„\ ,,J\ /J'',, {G = O(N)), (23b) 



Z r 



(2d 



l)- 1/2 P xy V xy Vj y pT (G = Sp(2iV)), (23c) 



where V xy is the Goldstone mode over Ai = U(2n)/ 
Sp(2n) or U(2n)/0(2n), resp., and P xy is the massive 
mode over Af/A4. Then we expand the Goldstone mode 
Uxy — exp(dx~y)U x up to quadratic order in the di- 
rectional derivatives d x - y , and Gaussian-integrate over 
the massive mode P xy = exp(X xy ) by retaining up to 
quadratic order in X xy . We reinstate the lattice spacing 
a and employ dimensionful continuum notations 



aM, £ d% -> W, Y, - ^" ' / 



Then the resulting effective partition function reads 

' N 



Z(M) 



VV exp 



M 



d a x 



HtVU^VU (24) 



NV2d l a i- dRe tr MU + Ca 2-d {tv V[/) - 



where U = V, V JV T , VV T for G = U(iV), 0(A), 
Sp(2A), respectively, and VV = Y[ x ^dV x . The coef- 
ficient C of the singlet (rf) kinetic term is a constant 
of O(l) depending on d and G B, but we suppress its 
explicit form as r/ should decouple from the Goldstone 
sector if the U(1)a anomaly were properly taken into 
account. Consequently the pion decay constant squared 
and the chiral condensate (the coefficients of the first and 
the second terms) are derived from the first principle. 

Conclusion. We have rigorously derived three classes 
of low-energy chiral Lagrangians from microscopic theo- 
ries: non-Abelian random flux models, or lattice gauge 
theories with staggered fermions in the strong coupling 
limit. These nonlinear a models have symmetric spaces 
M = U(ra),U(2n)/Sp(2n),U(2ra)/0(2n) as their target 
manifolds, depending on the complexity, reality, and 
pseudoreality of the defining representations of the gauge 
groups G = U(A),0(A),Sp(2A), respectively. These 
results are anticipated from Vafa-Witten argument |l2| ] 
or equivalently from the anti-unitary symmetries of the 
staggered Dirac operators in Eq.(Q), whose matrix ele- 
ments are trivially complex, real, and quaternion real 
for these three gauge groups. For the last two cases of 
the gauge groups, the anti-unitary symmetries of lattice 
Dirac operators are interchanged from those of the con- 
tinuum Dirac operators, due to the absence of the charge 



conjugation matrix acting on spinor indices, that squares 
to —1 14 1 . The puzzle of if and how the crossover be- 
tween different classes of effective theories could occur, 
originally posed for staggered fermions in the adjoint of 
SU(2) (= fundamental of 0(3)) and in the fundamental of 
SU(2) (m fundamental of Sp(2)), remains unsolved. Al- 
though the Wess-Zumino term in d = 2 dimensions is of 
the same order as the kinetic term, our models involving 
staggered fermions do not yield it. 

Verbaarschot |l5|] has conjectured that the three pat- 
terns of chiral symmetry breaking induce the spectral 
fluctuation of the Dirac operators to obey three univer- 
sality classes of chiral random matrices fllql . If our pa- 
rameters are restricted to the 'ergodic' domain where 
the size of the lattice is much smaller than the Comp- 
ton length of pseudo-Goldstone bosons, the chiral La- 
grangians are dominated by their zero modes, i.e. become 
finite-dimensional integrals involving only the mass terms 
in Eq.(|24|) [fL7| . As these three 'finite- volume' partition 
functions are known fl6| , |l8|| to be equivalent to chiral ran- 
dom matrix ensembles at (3 = 2, 1,4 in the large matrix- 
dimension limit, we have also proven his conjecture. 
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